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I. INTRODUCTION 



With the advent of the LHC, it is possible that in the not- 
too-distant future we will confirm supersymmetry (SUSY) as 
the next level of particle physics. Moreover, based on the pat- 
tern of superpartner masses, we may be able to conclude that 
SUSY breaking is of gauge-mediated type In this case 

it becomes probable that dark matter is due to gravitino LSPs 
iH]. Gravitinos can be generated in a variety of ways; via 
thermal scatterings |4-8], via inflaton and moduli decay 
ifTol [Till , and in the case of GMSB models, via decay of the 
energy density in the SUSY breaking hidden sector 1 12- 16]. 
In all cases a fundamental parameter is the reheating temper- 
ature Tr. Thermal gravitino dark matter generally requires a 
large Tr > 10^ GeV |17]. Production via inflaton decay, on 
the other hand, is strongly dependent on the branching ratio 
of inflaton decay to gravitinos lfT8l - l2lll . In the case of pro- 
duction via decay of a scalar field in the GMSB sector, the 
gravitino density will depend on both the initial amplitude of 
the scalar field oscillations and on Tr. 

It is generally not possible to obtain direct information 
about Tr from cosmological observations. However, it may 
perhaps be possible to deduce Tr from information obtained 
at the LHC. The point here is that Big-Bang Nucleosynthe- 
sis (BBN) can place strong constraints on the NLSP 1 17l l22ll . 
With enough information about NLSP couplings, it is possi- 
ble to calculate the thermal freeze-out density of NLSPs and 
to compare this with BBN constraints. If the thermal freeze- 
out density is too large, the simplest explanation in the context 
of the MSSM is that the thermal freeze-out NLSP density is 
diluted by a low reheating temperature below its freeze-out 
temperature. In this case we would need a mechanism to gen- 
erate the gravitino density with a low reheating temperature, 
Tr ^ 10 GeV. In the following we will consider whether such 
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a low reheating temperature can be consistent with gravitinos 
from the GMSB sector 

The SUSY breaking sector generally requires a gauge sin- 
glet field, in order to give masses to gauginos. At the end 
of inflation, its scalar component starts coherent oscillations 
around the minimum of its potential, due to the expansion 
of the Universe. A huge amount of energy is accordingly 
released, and typically a large spectrum of particles is pro- 
duced. This mechanism created much interest immediately 
after its discovery [23], as it can lead to cosmological difficul- 
ties which are related to either BBN [24] or to the overpro- 
duction of gravitinos lEsll . The latter comes not only from the 
thermal production through scatterings in the primordial bath, 
but in gauge mediation also from the decays of particles into 
gravitinos, which in this class of model are the LSPs. 

In this paper we investigate the possibility of dark matter 
gravitinos from the metastable SUSY breaking hidden sector, 
in particular for the case of a low reheating temperature ll26ll . 
In order to be able to explicitly calculate the gravitino density 
for all field amplitudes without cut-off, we focus on the case of 
a perturbative O'Raifeartaigh hidden sector In particular, we 
focus on models of gauge mediation with metastable vacua, 
which have been discussed by Muray ama and Nomura in [2^ . 

In this framework, gravitino production has been recently 
discussed by several authors (for instance, in lfT4ll . ifTsIl and 
^6]). It was found that the right amount of dark matter grav- 
itinos can be produced by the decay of the supersymmetry 
breaking field. Here we perform a study in the same direc- 
tion, though there are two important differences. First, in ifisll 
the minimum Sc of the scalar potential when the potential is 
dominated by Hubble corrections is independent of the Hub- 
ble scale and therefore constant in time. Once the S mass term 
dominates the potential, the field begins coherent oscillations 
with initial amplitude Sc. On the contrary, we show in Section 
im that the dynamics is quite different. The minimum of the 
potential when Hubble coiTections are important is time de- 
pendent and evolves towards the time-independent minimum. 
The S field tracks this minimum until Hubble coiTections are 
small and then enters coherent oscillations. This dramatically 
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alters the resulting gravitino density. Our main results are (i) 
the coherent oscillation amplitude and gravitino density are 
extremely sensitive to the parameters of the hidden sector and 
(ii) the reheating temperature can easily be as low as the BBN 
lower bound, Ir ~ 1 MeV, and still be consistent with grav- 
itino dark matter from the decay of the GMSB scalar. 

The inclusion of an O'Raifeartaigh hidden sector is the sec- 
ond feature of this work, since it is the first time that grav- 
itino production in GMSB with metastable vacua has been 
studied using a completely perturbative model rather than an 
effective model with a cutoff. In the latter case, the grav- 
itino abundance and Tr are dependent on this arbitrary cut- 
off (1 14-16]). By introducing an explicit hidden sector of the 
O'Raifeartaigh type, both A and the reheating temperature are 
explicitly determined as functions of the parameters of the the- 
ory. This means that the model at hand contains only quanti- 
ties, such as the coupling constants, which can be directly con- 
strained via cosmological observations or at colliders, since 
they determine the mass spectrum in the observable sector. 

The paper is organised as follows. In Section HIl we study 
the dynamics of the SUSY breaking field S and the ampU- 
tude of its coherent oscillations. In Section |III] we calculate 
the gravitino number density and the required reheating tem- 
perature under the assumption that gravitinos are the princi- 
pal constituents of dark matter In Section IV we present our 
conclusions. In the Appendix we calculate quantum and su- 
pergravity corrections to the scalar potential V{S) and derive 
bounds on the parameters of the model and on the gravitino 
mass range. 
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FIG. 1 : Time-dependent shift of the minimum of the potential due to 
the expansion of the Universe. 
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FIG. 2: The displacement from the minimum due to the time evolu- 
tion of Sq. 



II. COHERENT OSCILLATIONS OF THE SUSY 
BREAKING SCALAR FIELD IN AN EXPANDING 
BACKGROUND 

As already pointed out in the Introduction, our study of the 
dynamics of the 5-field differs from what is assumed in the lit- 
erature where the initial amplitude Sc of the scalar field 
oscillations is assumed to be determined by the Hubble cor- 
rections (with the Kahler potential normalized in units of the 
Planck mass). 



5kM ,^ 



(1) 



which imply \Sc\ = SkM/l^'K?. This assumes that the inflaton 
dominates the energy density of the Universe when the S field 
starts oscillating. We will also include this contribution in our 
analysis. However, in addition to ([U, the scalar potential con- 
tains terms which are not //-dependent. As we will see, this 
implies that the minimum 5'o(f ) is time-dependent in the pres- 
ence of Hubble corrections. The S field tracks the minimum 
as it decreases towards the time-independent (// = 0) mini- 
mum. This generates a displacement of S from the minimum, 
as shown in FiglT] As H decreases, either the quadratic or the 
linear term becomes dominated by the H — factor. So the 
minimum of V is //-dependent until both \S\^ and (5 + 5^) be- 
come //-independent. The displacement from the minimum at 



this time, 55, then determines the initial amplitude of the S os- 
cillations about the time-independent minimum ^o, as shown 
in FiglH 

This scenario assumes that the S field is not thermalized^ 
For T/{ less than the messenger mass M, the scattering of S 
particles by thermal MSSM particles will be due to messen- 
ger loops, with rate Fj « a^T^ /M'^, where a.; accounts for 
coupling and numerical factors. Comparing with the expan- 
sion rate at T < Tr, H ^ T^/Mpi, the hidden sector will not 
be thermalized as long as Tr < {M/asMpiY^^M. We will see 
later that this is easily satisfied in the models of interest to us 
here. 

We first discuss the dynamics of the GMSB sector gener- 
ally, then we will apply it to the specific case of the model 
defined in Section |III] Let us consider a potential of the gen- 
eral form 



V 



I 



mlS^-aS + 3H 



2 



cS 



(2) 



where the scalar field is assumed to take real values without 
loss of generality (with \S\ S/Vl for canonical normal- 
ization). When H decreases from a large value, either the 



The case of a thermalized hidden sector was discussed in |28 
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FIG. 3: Evolution of 55 in the potential with minimum at 652- 



linear or the quadratic term becomes dominated by the H- 
independent term. We define H\ and H2 by 



H = Hi 
H = H2 



when 
when 



a = 3cH\ 
ml = 2,H^. 



(3) 
(4) 



In the case H2 > Hi, the quadratic term becomes H- 
independent first. The potential when H2 > H > Hi is then 



V'. 



1 



mW - 3H^-cS. 



(5) 



On the other hand, if Hi > H2 then it is the linear term 
which becomes //-independent first. The potential when Hi > 
// > //2 is then 



H^S^ - aS. 



(10) 



The minimum is at 5o = a/3H'^. In this case Eq.O becomes 



55 + 55 = -y-3//255, 



which implies 



55i 



3fl 

2^2' 



(11) 



(12) 



The shifts (|9]l and (fT2l i give the initial amplitude of the os- 
cillations of S about the //-independent minimum in the two 
cases, where oscillations begin once H Hi and H Ri H2 re- 
spectively. 



III. DECAY OF THE SUSY BREAKING FIELD INTO 
GRAVITINOS AND CONSTRAINTS ON Tr 

A. A model of GMSB with metastable vacua and a hidden 
sector of the O'Raifeartaigh type 



The minimum is then 5'o(f ) — 3cH^/my The field S will track 
this time-dependent minimum. We define 5 = 5o + bS, where 
55 is a perturbation. Then the scalar field equation for S, 



S + 3HS- 



dV_ 
35' 



becomes^ 



...... dV 

55 + 3//55 = - (5o + 3//5o) - ^ 

o5 



Sa 



352 



(6) 



55. (7) 



So 



If the background is dominated by the inflaton, the Hubble 
parameter H depends on the scale factor a{t) us H °^ a^^l^. 
Thus the above equation takes the form 



55 + 3//55: 



21 cH'^ 



ml5S, 



(8) 



and the displacement 55 evolves as if it had a potential with 
minimum given by 



552 = - 
This is illustrated in Fig[3] 



27c//4 



(9) 



^ We have checked that the higher-order terms in the expansion of the poten- 
tial are negligible thi'oughout our analysis. 



We consider a superpotential 127 

W = -i/S + XSX^ + mXY + kSff + Mff, (13) 

and a non-minimal Kahler potential for the SUSY breaking 
field 5, 



K=\S\'-^-^2+0{^]. 



A4 



(14) 



The above equations define a model of gauge mediation with 
metastable vacua. In particular, there is a negative quartic term 
for 5 in the Kahler potential generated by radiative correc- 
tions, an accidental U{1)r symmetry which is broken in the 
messenger sector, and an explicit mass term for the messen- 
gers. As shown in Appendix A, the breaking of the additional 
symmetry induces a linear term in the Kahler potential, which 
makes the SUSY breaking field deviate from the origin. Such 
an effect, coupled to the expansion of the Universe, is respon- 
sible for the oscillations of the field 5 and accordingly, for its 
decay into gravitinos. 

The superpotential written in Eq.(fT3Tl corresponds to an 
O'Raifeartaigh hidden sector which contains, besides 5, the 
singlet fields X and Y. These have an approximately canonical 
Kahler potential (up to terms suppressed by the Planck mass). 
The parameters X and k can be taken real and positive without 
loss of generality, and the scales /j, m, and M have dimension 
of mass. The messenger sector contains Standard Model mul- 
tiplets belonging to 5 + 5* representations of SU(5), and the 
messengers / and / are assumed to have a canonical Kahler 
potential as well. Due to the second term in (fl4] i. a positive 
mass ms — /J^/A for 5 is generated around the origin. The 
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gravitino mass range allowed by the metastable vacua is in 
principle very large, 1 eV m3/2 ^10 GeV [IztI] . However, 
it is shown in Appendix B that it is possible to constrain this 
interval, by virtue of the explicit O'Raifeartaigh hidden sector 
and phenomenological constraints. 

As discussed in Appendix A and shown in iIztIi . radiative 
corrections generate A as a function of the mass scale m and 
couphng X in the hidden sector, 

A 2 ^Tt^m^ 

A=^. (15) 

Along the direction X = 0,Y = 0, Eqs.(fT3]l and (fT4] i localize 
the global SUSY minimum for the scalar potential at 



M 



(16) 



Moreover, to avoid tachyonic messengers, we require that the 
condition 



> k// 



(17) 



is satisfied. In this case, we find a local minimum at 5 = / = 
/ = 0, which is a metastable vacuum. At this point in field 
space supersymmetry is broken, with — /j^. As we will see 
shortly, 5 is a flat direction at the classical level, but it is lifted 
around the origin by the radiative corrections to the Kahler 
potential that are generated by the loops of the fermions and 
scalars in the hidden sector and of the messengers. The mass 
matrix in the basis of the X and Y fermions and Vf/y is 

2XS m 



Ml 



m 



(18) 



while the F-terms of the chiral superfields, in the basis 
{X ,X\Y,Y^), produce the following scalar mass matrix ll27ll 



/ot2+4A.2|5|2 -2V 2A.m5^ \ 



Ml 



V 



-2V 
2XmS 




nr+4A,2|5p 




m 




2XmS 


2 

m 



(19) 



For the messenger sector, it can be shown that the scalar and 
fermion mass matrices M^' and MJjJ, in the bases (/,/^,/,/^) 
and can be respectively written as 

m;^ 

( \M + kS\^ -ki? + kXx'^^\ 

\M + kS\^ -kp^ + kXX^ 

-kf? + kXX^^ \M + kS\^ 

\^kfi^ + kXX^ \M + kS\^- ) 

(20) 



and 



M + kS 
V~\M + /t5 



(21) 



The radiative corrections generated by these particles shift the 
minimum of the scalar potential only slightly from the ori- 
gin. However, the shift due to these terms is important from 
a cosmological point of view, since it is responsible for the 
production of gravitinos through the decay of the S field, as it 
is shown in the following. 



B. Gravitino dark matter and reheating temperature 

The results of Section|II]and the GMSB model summarised 
above can be now applied to gravitino production. Here we 
calculate the gravitino number density as a function of the re- 
heating temperature Tr and the masses and couplings of the 
O'Raifeartaigh sector. The scalar potential of the field S is 
obtained from the superpotential ( fTST i and from the Kahler po- 
tential ( fT4l i. as a sum of a classical part and the corrections 
generated by quantum loops and by supergravity. Each of 
these contributions is discussed in detail in the Appendix. In 
particular, the one-loop corrections ll29ll to the Kahler poten- 
tial can be calculated by using the following expression POtl : 




2Mln 



M 
A2 



(22) 



which is generally valid for a renormalizable N — 1 theory. 
Here M| and My are the mass matrices of the chiral and vector 
superfield sectors respectively, with 

iMl)ij ^Wf^^'^'^Wkj or pik=mik + Xijk<t>\ (23) 



for the chiral sector and 



Mab — 2 ( ^A^Bi + X'gXAi ) , 



(24) 



for the vector sector (which is however absent in the model 
discussed here). By summing equations ( IA9I ) and ( IA26b . we 
obtain 



ViS) = VsusY- 
2A.V'* 



5k 
167i2 



3H 



M{S + S'') -3H 



f Vsugnl ' 

48n^M^ 



\S\' + 

2ky 



(25) 



where we have kept only the linear and quadratic terms in 5, 
in the direction corresponding lo X — Y — f = f — 0. These 
terms are dominant for sufficiently small values of the field. 

The scalar field can be treated as a modulus, with interac- 
tions given in ifTill . The decay into gravitinos is not helicity 
suppressed, and gravitino production from S decay dominates 
production from inflaton decay^. 

Dominant decay into a pair of gravitinos {S — > "^i^s^v) occurs 
if the cutoff A satisfies ifTsIl 



A< 8 X 10' 



GeV 



(26) 



' The mass scales here considered are consistent with the discussion in 



Ref.yJ. 
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The inclusion of the O'Raifeartaigh sector makes it possible 
to relate this directly to the mass scale m in the hidden sector 
using Eq.lfTSll. 



,14^2 



X"- GeV X 



m3 
1 TeV 



)"'(t 



GeV 



(27) 

This corresponds to an upper bound in the range 10** GeV- 
10^^ GeV for a gravitino mass between 1 keV and 10 GeV 
(the lower bound on m3/2 is discussed in Appendix B). 

The messenger mass scale M is related to the gluino mass 
m3 and to the SUSY breaking scale /j by ifTsIl 



M = 



47im3 



(28) 



Therefore, by using the relation between and 

= V3my2Mpi , (29) 

with the reduced Planck mass M/./ = 2.43 xlO^** GeV, we find 

'1 TeV\ / ^3/2 



Mw3.3 X lO'^A: 



\0.lJ 



mi 



J V 1 GeV 



GeV . (30) 



This impUes 10^ GeV < M < lO"* GeV for the above range of 
gravitino mass. For this range of messenger mass, the corre- 
sponding upper bound on the reheating temperature for which 
the S field is not thermalized is 1 TeV to 10^" GeV (using 

1/3 IS 

a/ - 1). 

Eq.(l30ll. together with dZTl l. gives the following approxima- 
tion for the coefficient of the linear term in the scalar potential 



2^/4 2ytV^ 



M2 



2kY 

M2 



(31) 



over the entire gravitino mass range considered. Similarly, for 
the coefficient of the quadratic term. 



5kY \ 



■ + 



2XV 
3n^m^ 



(32) 



One should also consider what happens when adding a con- 
stant term ^ m2/2Mp[ to the superpotential of the model. This 
is needed to tune the vacuum energy to zero in supergravity 
Recalling Eq.(fT3]l. we can write il4ll 



W 



-IJ^S + my2Mli + . 



(33) 



This generates both linear and quadratic terms in S through 
the factor jWp in the supergravity corrections ll3lll . 



dw 



W 



Mil ^^^J 



W dK 



-1 



-3 



Ml, 



(34) 



The authors acknowledge Kazunori Kohri for this idea. 



In fact, the contribution of /2^p; is 111 



Vsugra = &Xp{K/Ml,) 



.-3- 



\W\ 



M- 



1 + 



PI 

Mp,/ 

5k M 



Ml, 



167l2 Ml, 



6p^my2{S + S^) 



Mpi 



\5k Mfi'* 
4V3n^W' 



PI 



Eqs.(l3Tl) and (|32] | change accordingly: 
^ 2fcV 32\/37tV 



3H^-- 



3H^ 



1 



M2 



2XV 



3Ti^m^ 



5k MMpi 

15;t/ M 
W3n^Ml, 



(35) 



(36) 



(37) 



For the entire gravitino mass range here considered, the con- 
tribution of the additional term in dJTl l is always negligible. 
On the other hand, if A: 0. 1 the new term in ( [36] l dominates 
for certain values of 12, thus we take it into account as well. 
The scalar potential of interest is therefore the following. 



3H^ 



2XV 
3n^m^ 



+ 



5kM S 



-3//2 



2fcV 32\/37i:V 



5k 



MMpi 



(38) 



where it is assumed that S is real for simplicity, which will be 
the case on minimizing the potential if all parameters are con- 
sidered real. Equation ( l38b can thus be rewritten as follows. 



c2 , f 5ky 



\4V2%^M MpiJ 



-3H^ 



52 



5kM 



2 8V27i2 

Comparison of this expression with the potential 



V = ^m^S^ - fl5 + 3//2 f'y - c5 



(39) 



(40) 



gives the parameters which determine the dynamics of the 
field S, as discussed in SectionHJ 

We next review the calculation of the gravitino abundance 
from decay of a scalar field S. Oscillations of S start when the 
Universe is still dominated by the inflaton, thus T > Tr. Let 
nsosc be the number density of the scalar and Oosc the scale 
factor at that time. The number density ns{T) at a generic 
temperature is given by 



nsiT) 



a{T) 
H{Tr) 



nsoi 
2 



a{TR) 
a{TR)J \a{T) 



g{T)T' 



Hose J g{TR)T, 



3 osc 1 



(41) 
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since H °^ a^^l^ during inflaton domination. Recalling the 
expression for the Hubble scale during radiation domination, 
the above becomes 



ns{T)- 



%^ g{T)T^TRnsos 
90 Mil Hi, 



Therefore the number density to entropy density of S is 



^ 4M2, Hi, 



(42) 



(43) 



From the scalar potential Eq.([39]l. ps = which im- 

plies nsosc = Ps/tns = /2)'«S'5ast, where Sosc is the value of 
the field at the beginning of the oscillations. Therefore 



osc 



s 8Mj,Hl, 



(44) 



Since each S scalar decays into a pair of gravitinos, the grav- 
itino number density to entropy ratio is given by 



and 



Tr 



-4 



1 TeV 



'«3/2 



1 GeV 



1014 GeV/ V 10*^ GeV 



127i:a3 / ^3/2 



(52) 



By demanding that the gravitino is the principal constituent of 
dark matter, 



03/2/^2 «0 (0.1), 



(53) 



we then obtain the reheating temperature 7^ as a function of 
the parameters of the model 



;7i3 



A.14 /a3\4 

w0.6GeVx — i- ( — ) , 

kf^ \0AJ VI TeV 



VI GeV/ V1014 



GeV 



-7 



1 



127ra3 / 
5k^ \ niT, 



-3 



(54) 



"3/ 2 _ msTRSi 
s 



4Mj,Hl, 



(45) In the case Hi > H2, Hi, ss //j = m^/2i and Sosc = ^Si . We 
then obtain 



where the initial oscillation amplitude around the minimum, 
Sosc, is given by 5Si or 552 at the time when all H dependence 
becomes negligible. 

Consider first the case where H2 > Hi, namely where 
-a /3c and Sosc = 5^2. Therefore 



f^osc ~ f^i 



ns msTr bS2 
s 



mj, Hf 



(46) 



The parameters a and c and the mass ms are fixed by the mes- 
senger and O'Raifeartaigh hidden sectors. 



2X4/ 

4V2n^M 
5kM 



16\/37i:2 M " 
5k^ 



Using these we obtain 



ns 202571'^ 0m\m''TR 



I6V2 A-i^a^m^/jM^, 



where we have used Eq. 
/J. Therefore 



(47) 
(48) 
(49) 

,,21 , , (50) 

Dk-^ V m3 

and Eq. ( |29] l to eliminate M and 



127ra3 / »J3/2 



m3 



^«0.11xilf^)"\ 
s A.14V0.I/ VI TeV 



»^3/2 

1 GeV 



1014 GeV / V 10^ GeV 



Tr 



127ta3 f 'W3/2 \ 
5/t2 m3 y 



(51) 



«s _ 1822571^ k'^mlm^TR 
s 



256V2 ^^i^a^m^^j^l,, 



127ia3 f mi 2 
5k^ \ m$ 



(55) 



Therefore 

s 



3.9 X 10-3^ 



X 



1014 GeV 



Xio VO.I 

5 / Tr 



m3 



'/ VI 



^3/2 



and 



3/2 



/i2 w5.9x 10^ 



10^ GeV 



^ / a3 
?Jo VoT 

5 / Tr 



1 TeV / VI GeV 
^ 127ia3 / »J3/2 
5k^ \ niT, 



(56) 



/ "^3/2 y 

VI GeV/ 



1014 GeV/ VlO^GeV 



m3 
1 TeV 



127ta3 / OT3/2 
5k^ \ m3 



n 2 



(57) 



The reheating temperature Tr required for the correct density 
of gravitino dark matter is then 



Tr w 17 GeV X 



a3 



fc4 VO.I 



( "^3/2 \2 / 

VI GeV/ V: 



1014 GeV 



-5 



1 TeV 



127i:a3 / 'W3/2 
5k^ \ »J3 



(58) 



Eqs.(l54l) and ( |58] | are the main results of this article. We see 
that the model can accommodate a very wide range of Tr and 
still be consistent with gravitino dark matter from the GMSB 
sector. The most striking feature of these results is their ex- 
treme sensitivity to the parameters of the model, in particu- 
lar the O'Raifeartaigh sector coupling X. This means that the 
gravitino density from decay of the SUSY breaking scalar in 
the GMSB sector can account for dark matter for essentially 
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any value of the reheating temperature above the BBN bound, 
Tr > 1 MeV. 

Comparing our results with the previous results of ifTsll . 
which assumed that oscillations about the minimum began at 
Hoic ~ tns with S^sc ~ V^c (for real S), we find; 



»3/2 _ »3/2 
»3/2 _ »3/2 





X 






1 










^' 


cm| 




X 






1 


fl 


£7 






4 ' 





where n3/2/'s|£j- is the value given in 11151 



«3/2 
S 

X 



:2.,0-».*(^)^(- 



TeV 



7> 



»^3/2 \ 

1 GeV/ V 10** GeV / V lO^^ GeV 



< 1 , (59) 
> 1 , (60) 

-2 

(61) 



If Hi > H2, \a/cm^\ > 1 and so there is a strong enhancement 
of the gravitino abundance relative to the amplitude obtained 
in Us]. Similarly, if H2 > Hi then there is a strong suppres- 
sion of the gravitino abundance relative to lITsIl . 

Note that while in the existing literature Tr is necessarily 
written in terms of an arbitrary cutoff ifTsll . in our analysis 
there is a fully defined perturbative hidden sector, eliminating 
the cut-off in favour of the masses and couplings of the model. 
Therefore what could have been a limitation, namely focusing 
on a specific O'Raifeartaigh -type sector, is in fact an advan- 
tage. Since there is no cut-off, there is no need to constrain the 
parameters of the model in order to make it consistent with the 
validity of the effective theory. 



IV. CONCLUSIONS AND OUTLOOK 



superpotential coupling of the O'Raifeartaigh sector As a re- 
sult, it is easy to account for gravitino dark matter with an 
arbitrarily low reheating temperature. (There are also con- 
straints on 111^/2 ^nd X from gravitino free-streaming, but as 
shown in Appendix C these are easily satisfied.) 

This could be significant for the cosmology of GMSB mod- 
els. It is possible that the NLSP could be discovered at the 
LHC and its properties established. In particular, the thermal 
freeze-out density and decay rate of the NLSP may turn out 
to be inconsistent with BBN. Should the LHC discover a cos- 
mologically problematical NLSP, a low reheating temperature 
would be the simplest way to solve the problem, by suppress- 
ing the NLSP density via entropy release prior to Tr. This 
requires that Tr is sufficiently small compared with the freeze- 
out temperature of the NLSP, Tr < Mnlsp /20 < 50 GeV. In 
this case the GMSB sector would be a prime candidate for the 
origin of the gravitino dark matter density in the presence of a 
low Tr. (However, the reheating temperature required to pro- 
duce the correct density of gravitino dark matter is very sensi- 
tive to the parameters of the model, so obtaining a low enough 
reheating temperature to dilute the NLSP density while not af- 
fecting BBN may require a degree of coincidence.) To com- 
plete the model a source of baryogenesis consistent with a low 
Tr would be required. A natural possibility would be Affleck- 
Dine baryogenesis. Indeed, for m^i2^ 1 GeV, Affleck-Dine 
baryogenesis combined with Q-ball decay in GMSB could si- 
multaneously account for both gravitino dark matter and the 
baryon asymmetry ll33ll . but for smaller m3/2 it would provide 
only the baryon asymmetry. 

As an additional remark, we note that the 
O'Raifeartaigh GMSB sector is a well-defined perturba- 
tive model. Therefore the decay of the SUSY-breaking scalar 
field is addressed without any cutoff, nor with any stringent 
constraints on the mass scales nor on the coupling constants. 



In this paper we have considered the dynamics of the SUSY 
breaking scalar 5 in a GMSB scenario with metastable vacua 
and the production of gravitino dark matter through its decay. 
Our results for the cosmological evolution of the scalar field 
and the resulting gravitino density are significantly different 
from previous investigations. 

We have shown that since the Universe is expanding, the 
minimum of the potential Sq is time-dependent. The S field 
tracks the minimum, which generates a displacement hS. 
Once the minimum becomes H independent, S begins coher- 
ent oscillations about the minimum with initial amplitude de- 
termined by the displacement. This produces a very different 
S oscillation amplitude and so gravitino density as compared 
with previous analyses [15]. By considering a generic poten- 
tial V{S), we have shown that there are two possible values 
of 55, depending on whether the quadratic or the linear term 
first becomes //-independent. The resulting gravitino density 
can be highly suppressed or enhanced as compared with the 
previous estimate of fTs'l, depending on the parameters of the 
model. A striking feature of the gravitino density is its ex- 
treme sensitivity to the parameters of the model, with the re- 
heating temperature having a A-'" or A,^^ dependence on the 
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Appendix A: Radiative corrections at one loop and effective 
potentials 

In this section we calculate the radiative corrections to the 
scalar potential V{S) and to the Kahler potential of S. The 
supergravity contribution is obtained accordingly. 
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1. Contribution of global SUSY 

The potential Vsusy includes the classical contribution Vq 
and the effective potential Vg/f that is generated by loops of 
X, Y and of the messengers / and /, 



Vsusy (S) = Vo{S) +Veff{S) ^ V/;{S) +V(;'{S) 
+V'{S) + V!:US). 



(Al) 



The classical scalar potential of the hidden sector, Vq has the 
following form 



+2Xm{SXY'^ +S<X''Y)+4X^\X\^\S\^ , 



\Y\') 
(A2) 



which is given by the F-terms derived from Eq.(fT3b. 

In general, the scalar and fermion loops correspond to the 
Coleman-Weinberg potentials i29ll 



Tr 



for the scalars, and 



y. 



eff- 



6471:2 



Tr 



(|M,|^)^ln^ 



(A3) 



(A4) 



for the fermions. The minus sign comes as usual from the 
Fermi-Dirac statistics. The mass matrices My and are 
given by equations ( fTsT i and (fT9] l respectively, and R is an ul- 
traviolet cutoff. 

The formulas (|A3t and (IA4b give the following Coleman- 
Weinberg potential for the hidden sector. 



v^ff{s) = v'^{s)+2vi;{s)'^ 



\S\' 



107l2m4 ' 



(A5) 



where we have dropped a constant and kept the dominant 
terms in XjA^/m'^ and in |5|2. The above equation agrees 
with [27]. The positive mass term, which for small values 
of S dominates over the quartic term, confirms what has been 
claimed in Section |III] By adding now Eq. (|A5l l to ( IA2b . the 
scalar potential of the hidden sector can be recast as follows, 

^^iUX^\x\^-X^l\x^+x^^)+m\\x\'' + \Y\^) + 

+2Xm(5XF^"+5+Xl'y)+4A.2|X|2|5|2 + 



37i2m2 



107l2m4 



(A6) 



Considering now the messenger sector, one can act in com- 
plete analogy. The F-terms of the messengers give the fol- 
lowing tree-level potential: 

y^' = k{-^? + XX''-) ff + k{-^? + Xx^)f-f + fc'l/p l/T . 

(A7) 

The mass matrices of the messenger scalars and fermions are 
given by equations (|20] | and (ISTT i in Section |III] By summing 



the two contributions as before, and keeping only the dom- 
inant terms in kjj^ /fvfi, in X^P' /m'^ and in S, we obtain the 
Coleman-Weinberg potential of the messengers. 



3,, 2 



5ky 

' 327t2M2 



5A:> 
167i2 



M 



X 



t2^ 



(5 + 5'^ 



(52 



-5^2) 



247t2M6 



X2 



^'X\xY)+X^?{x^+x'^) 



(AS) 



where we have dropped an unimportant constant. The factor 
of 5 comes from representations of SU(5). This expression 
agrees with the literature lHH^l for X = F = 0. Notice that 
the messenger loops generate a mass term for the field S as 
well. 

By adding ( IA7b and (IA8b . the scalar potential of the hidden 
and messenger sector along the direction X = Y= f= f = 
in the field space, which is our case of interest, can be written 

as 



3, ,4 



Vsusy (S) ■ 



|5|2- 



4, ,4 



327i2m2 



(52 + 5^2). 



487i:2m6 y ' ' 327i:2m7 

Il07l2„j4 + 647I2MV' ' 



(A9) 



This will be added to the supergravity corrections of the next 
section, in order to obtain the full potential of the model. 



2. Kahler and supergravity corrections to the potential 

Let us first focus on the Kahler potential. In the case at 
hand, Eq.(l23b is rewritten as = (1 — hik){M + kS) for the 
messengers. Therefore Eq. (l22l l can be rewritten as follows. 



Am — — 



327l2 



Tr 



-m|2 



|m;«|2 In^ 



- 1 



(AlO) 



where, as before, the factor of 5 corresponds to SU(5) repre- 
sentations. At the lowest order in kpP jM^ and for small values 
of 5, we obtain 



--(- 

2 \M 



5M2 
16%^ 



|5|2(5 + 5+)- 



1 

6 \M 



|5|2(52 + 5+2) 



(All) 

as it was found also in fT3\. On the other hand, loops of the 
superfields X and Y in the hidden sector generate the following 
Kahler potential, expanded in S and in X/j^/m^, 



(A12) 
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The above agrees with the statement reported in tZm . namely 
these corrections correspond to the Kahler potential 



with 



8jt2 



. 7 3k in 

A = 



?l4 



(A14) 



since X^/Sn^ <C 1 even for a strong coupling ^ « o(l). 

To embed the model in supergravity, we can use the general 
formula ifsl 



Vsugra = exp{K/Mji) 



dW W dK\ /dW W dK 



Mil 



with the Kahler derivative of the superpotential W 
'dW W dK 



DW = 



d(^a Ml, d<^a 



(A15) 



(A16) 



At this point, it might be instructive to show the effect of the 
quartic term in S in the Kahler potential at the tree level, be- 
fore including the radiative corrections from the loops of the 
messengers and of the hidden sector. Eq.( IA15l ) for the SUSY 
breaking field is rewritten in our case as 



Vsugra{S)^IJ^exp{K/Mli) 



1 - 



\sl 

A2 



X 1+2 



^ \S\' 



M 



PI 



PI 



\S\' 



M 



PI 



(A17) 



with the Kahler potential ( IA13I ). At the lowest order in \S\^, 
this is expanded as 



|5p 



7 151' 



A2 4 A2m2, 



1 



X^ |5|2 7A.4 \S\* 



37t2 



(A18) 



The term proportional to I^P has a positive coefficient, con- 
sistently with the literature llisll . 

If the radiative corrections ( lAl 11 1 and ( IA12I ) are now taken 
into account, the dominant terms are 



k'^'\s) = i6^\s)+k^,^\s)^\s\^ 

5k' |5|2 



327t2 M 



5k 

167l2 

A.4 



M{S + S^)- 



127t2 m2 



(A19) 



thus the following dominant terms are generated. 



Vsu 



TTT ' 



15V- 



M 



PI 



64714 



M 
Mpi 



"37i2 \ m J 



5V 



7(5 + 5^) 2^ 

167t2^ '\M m2 



M 



PI 



(A20) 



where we write only the terms which are linear and quadratic 
in S. This agrees with dAlSb in the limit — > 0. We can further 
assume that 



25k^ 
647? 



M Y 



37t2 



m J 



(A21) 



therefore the dominant supergravity corrections can be rewrit- 
ten as 



X4 



5k 



167:2-(^ + ^nM-^M2, 



M 



(A22) 

Clearly, the sign of the term that is linear in S is now deter- 
mined by the coupling constant k. In particular, a form of the 
potential such as 



Sugru \ 



+c^\S\^~d^{S + S^) + ... 



(A23) 



would be recovered if the following condition is satisfied: 



^< v2— 
Mpi 



my2 ^ 5.77 



m3 



a3 



(A24) 



The above equation directly relates the gravitino mass mo, /2 to 
the gluino mass ms and to the strong coupling constant as. It 
is calculated by using Eq.([30l) in Section|IIIl However, ( IA24I ) 
is not consistent with gauge mediation, where the gravitino is 
the LSP. In particular, the framework that is used in this paper 
allows gravitino masses 10 GeV, thus we conclude that 

also the coefficient of the linear term in (IA22b is identically 
positive. 

The supergravity corrections to the scalar potential of the 
model can now be written by adding ( IA22b to the potential 
for S which is generated by the inflaton, whose energy density 
dominates the Universe before the field starts oscillating ifisll . 



y(5)«e^(3//2M^)«3//2 



5k 
T67t2 



M(5 + 5^) 



(A25) 

(where K has been normalised in units of Planck mass). This 
provides, at the lowest order in the SUSY breaking field. 



VsugraiS) 

X' 



W3//2 

|5|2-- 



|5|2 
5k 

167t2 



5k 

167? 



M{S + S'' 
M 



'Mil 



(A26) 



namely with the supergravity corrections which, together with 
Eq. (|A9l l, are used in Section |IIT] to calculate the gravitino 
abundance and the reheating temperature. 
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Appendix B: Consistency of the Coleman Weinberg potentials 
and constraints on the gravitino mass range 

In this appendix we first note that the coupling constants of 
our model X and k and the mass scale m are constrained by 
requiring consistency of the series expansion of the Coleman- 
Weinberg potential. The scalar potential Eq.dZSl) is obtained 
by expanding in terms of Xfj^/in^, so the condition 



^<1, 



must be satisfied. Therefore 



m > 



(Bl) 



(B2) 



Next we derive further constraints on the mass scale m of the 
hidden sector, and show that the gravitino mass range which 
is consistent with the model can be accordingly constrained. 
First, using Eq. ( |29] |. the lower bound (IB21 l can be recast as 



m > = 2 X 10\'/2 GeV 



'W3/2 
1 GeV 



(B3) 



In Section Un] we assume that the gravitino is the main decay 
product of the field S when this rapidly oscillates around the 
minimum. Accordingly, the upper limit Eq.(l27b holds: 



If GeV X 



1 TeV 



-1 



VI GeV/ 



(B4) 

Clearly, ( |B3l l and ( IB4I ) must be simultaneously satisfied. This 
means that the formulas ( |54] i and (|58] l are valid only for cer- 
tain gravitino mass ranges, depending on the values of X. In 
particular, 



1: 



10-2 = 

no-1 



1 MeV< m3/2< 10 GeV, 
> 1 keV<OT3/2< 10 GeV, 



(B5) 
(B6) 



where we have assumed only that ;«3 sa 1 TeV. Remark- 
ably, since these results follow from the consistency of the 
Coleman-Weinberg potentials and from phenomenology, they 
do not rely on any approximations. 

Let us finally show why the coupling k can be consistently 
of order unity. In analogy with Eq.( IB2b . the effective potential 
I for the messenger sector is valid if 



ki? , 



(B7) 



By recalling Eq.(l30b, i.e. 



Mw3.3x 10'^ GeV X A; 



1 TeV 



»^3/2 

1 GeV 



(B8) 



we can substitute the above equation and ( |29] l in ( |B7| i to finally 
obtain 



10- 



Ik 



0.1 \- / ma y / m3/2 ^ 
a3 / VI TeV/ VI GeV/ 



(B9) 



This condition is satisfied for k ^ 0.01 in the entire gravitino 
mass range in equations iB5i and ( IB6I ). 

Appendix C: Constraints from gravitino free-streaming 

One should also take into account the free streaming length 
XfS, namely the distance that the gravitinos produced through 
the decay can travel until they become non-relativistic. If we 
want them to be Cold Dark Matter (CDM), this quantity must 
be smaller than (lOOkpc). 

The free-streaming length for gravitinos from S decay was 
derived in ifTsll . 



V.^lkpcJ— 



3/2 



1014 GeV / 



(CI) 



where g^, is the number of relativistic degrees of freedom at 
the time of decay. Using A = \/3nm/X^, this becomes 



1 kpc 



^3/2 

100 V IGeV 



lO'^GeV 



(C2) 

Since the free streaming length must be smaller than 100 kpc, 
one obtains a lower limit on the gravitino mass, 

-^/-«-2^4S)HtoT^)°^^- 

This depends on m and X. Using the lower bound on m, 
Eq. (IB3l l, we obtain a lower bound on the gravitino mass from 
the requirement that gravitinos behave like cold dark matter 



m3/2>2.6xlO-eVx^(^)' 



(C4) 



From this we see that free-streaming imposes only a weak 
constraint on the gravitino mass, unless X^l. 
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